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In a ferromagnet, an anomalous-Hall heat current, given by the off-diagonal Peltier term a^y, 
accompanies the anomalous Hall current. By combining Nernst, thermopower and Hall experiments, 
we have measured how a^y varies with hole density and lifetime t in CuCr2Se4_xBrx. At low 
temperatures T, we find that a^y is independent of r, consistent with anomalous-velocity theories. 
Its magnitude is fixed by a microscopic geometric area A ^ 34 A^. Our results are incompatible 
with some models of the Nernst effect in ferromagnets. 
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In a ferromagnet, the anomalous Hall effect (AHE) is 
the appearance of a spontaneous Hall current flowing par- 
allel to E X M where E is the electric field and M the 
magnetization Karplus and Luttinger (KL) 0] pro- 
posed that the AHE current originates from an anoma- 
lous velocity term which is non-vanishing in a ferromag- 
net. The topological nature of the KL theory has been 
of considerable interest recently H, Hi IE 111 - Experimen- 
tally, strong evidence for the dissipationless nature of the 
AHE current has been obtained in the spinel ferromagnet 
CuCr2Se4_xBrx. Lee et al. reported that, despite a 
1000-fold increase in the resistivity p induced by varying 
the Br content x, the anomalous Hall conductivity (nor- 
malized per carrier and measured at 5 K) stays at the 
same value, in agreement with KL's prediction. A test 
of the anomalous velocity theory against the AHE in Fe 
has also been reported 7]. 

It has long been known that an anomalous heat cur- 
rent density J*^ also accompanies the AHE current in the 
absence of any temperature gradient In principle, 

can provide further information on the origin of the 
AHE, but almost nothing is known about its properties. 
A weak heat current is a challenge to measure. Instead, 
one often performs the 'reciprocal' Nernst experiment in 
which a temperature gradient — VT produces a trans- 
verse charge current, which is detected as a Nernst elec- 
tric field Eat parallel to M x (— VT). However, in previ- 
ous Nernst experiments on ferromagnets |E E, 13 ) 
was not found because other transport quantities were 
not measured. Combining the Nernst signal with the 
AHE resistivity p'^y and the thermopower, we have de- 
termined how the transport quantity axy relevant to 
varies in CuCr2Se4_xBrx as the hole density Uh and car- 
rier lifetime r are greatly changed under doping. We 
show that axy has a strikingly simple form, with its mag- 
nitude scaled by a microscopic geometric area A. 

We apply a gradient — VT| |x to an electrically-isolated 
sample in a magnetic field H||z. Along x, the charge cur- 
rent driven by — VT is balanced by a backflow current 
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produced by a large which is detected as the ther- 
mopower S = Ex/\VT\. Along the transverse direction 
y, however, both and — VT generate Hall-type cur- 
rents. In general, the charge current in the presence of 
E and - VT is J = cr-E -I- a ■ {-VT) with cr and a the 
electrical and thermoelectric ('Peltier') conductivity ten- 
sors, respectively. Setting Jy = 0, we obtain the Nernst 
signal Cat = Ey/\VT\ = pa^y + PxyOt where a = Uxx IS- 
Hence, as noted, the Nernst signal results from the 2 
distinct y-axis currents ayx{—VT) and GyxEx. In a fer- 
romagnet, the former is our desired gradient-driven cur- 
rent, whereas the latter comprises the 'dissipationless' 
AHE current and the weak ordinary Hall current. 

In terms of the thermopower S — pa and Hall angle 
tan^^f — Pyx/ p, we may express axy as 

paxy = eN + Stand H- (1) 

Hence, to find axy, we need to measure cn, S, pxy and p. 
Knowing axy, we readily find the transverse heat current 
= OiyxEx since ctyx = ayxT by Onsager reciprocity. 

The spinel CuCr2Se4 is a conducting ferromagnet with 
a Curie temperature Tc ^ 450 K. Because the exchange 
between local moments in Cr is mediated by superex- 
change through 90° Cr-Se-Cr bonds rather than the car- 
riers, Tc is not significantly reduced even when the hole 
population rih drops by a factor of 30 under Br doping 
{M at 5 K actually increases by 20%) H Using 
iodine vapor transport, we have grown crystals with x 
from 0.0 to 1.0. As a; increases from to 1, the value of 
p at 5 K increases by ~ 10'^, while Pxy/uh increases by 
~ 10^ 0. The tunability of Uh and the robustness of M 
under doping make this system attractive for studying 
charge transport in a lattice with broken time-reversal 
symmetry. The behavior of p, M, and p'^y vs. x are 
described in Ref. @. 

Figure n shows profiles of ca? vs. H at selected T in 2 
samples with x = 0.1 and 0.85 and Tc = 400 and 275 K, 
respectively. As noted above, eN{T, H) is the sum of two 
terms both of which scale as M. The magnitude |eAr| ini- 
tially increases as H rotates domains into alignment and 
then saturates to a constant for H > H^, the saturation 
field. The sign of cat - negative in all samples - refiects 
the sign of the dominant term . 
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FIG. 1: Curves of the measured ejv = Ey/\\/T\ vs. in 
CuCr2Se4-xBrx with x — 0.1 (left panel) and 0.85 (right). In 
the ferromagnetic state below Tc, ejv saturates to a constant 
when H exceeds Hs reflecting the M-H curve. The scaling 
factor Qs increases rapidly as T increases from 10 K to Tc- 
In the right panel, ejv continues to scale as the M-H curve in 
the paramagnetic regime (275-400 K). 



In the sample with x = 0.85, the curves above Tc show 
that the scaling also holds in the paramagnetic regime 
where the susceptibility has the Curie- Weiss form x ~ 
1 / (T—Tc) in weak H . In analogy with the Hall resistivity 
Pxy — Rq^qH -\- Rs^qM , with i?o and Rs the ordinary and 
anomalous Hall coefficients, respectively, it is customary 
to express the scaling between the e^q-H and M-H curves 
by writing 

eN = QoMofl" + QsMoM. (2) 

For T < Tc in all samples, the Qq term cannot be re- 
solved, so that ejy ~ Q^fioM. Moreover, below 50 K, 
M changes only weakly with x (by 20% over the whole 
doping range) , so that the saturated value of the Nernst 
signal e^* differs from Qs by a factor that is only weakly 
X dependent. 

The Nernst signal has very different characteristic be- 
haviors below and above Tc- As an example. Fig. [21 
shows e^* measured at 2 T in the sample with a; = 1.0 
(Tc = 210 K). Between 5 and 100 K, e^* increases lin- 
early with T- Above 100 K, e^* rises more steeply to 
a sharp peak 200 K, and then falls steeply above Tc- 
As noted, in the paramagnetic regime, the Nernst signal 
matches the behavior of M as a function of both T and 
H- Figure El shows that the T dependence of bat closely 
follows that of M — xH (both are measured at 2 T). 
The experiment shows that, in a gradient, fluctuations 
of the paramagnetic magnetization lead to a significant 
transverse electrical current that is proportional to the 
average magnetization (this has not been noted before. 



FIG. 2: The T dependence of the Nernst signal cn (solid 
triangles) measured at 2 T in the sample with x = 1.0. Above 
Tc, ejv is compared with the paramagnetic magnetization M 
at 2 T (open circles). 

to our knowledge). We express the proportionality as 

a^y^l3M {T>Tc), (3) 

where P is only weakly T dependent (it decreases by 5% 
between 250 and 400 K). The parameter /? plays the im- 
portant role of relating the magnitudes of the paramag- 
netic M and the transverse electronic current (through 
the Nernst signal). Its very small value (/? ~ 2 x 10^^ 
at 250 K) reflects the strikingly weak coupling be- 
tween the fluctuating M and bat in ferromagnets, com- 
pared with vortex flow in a superconductor [T2 . [l5l |. With 
the growth of long-range magnetic order below Tc, Eq. 
ceases to be valid. 

In the ferromagnetic state, we restrict our attention to 
the regime below 100 K where e|^' is nominally linear 
in T- Figure shows curves in this regime for the 5 
samples studied. The slopes of the low-T curves are not 
monotonic in x. As a; is increased from 0.1, the slope at- 
tains a maximum value at x = 0.25, and then decreases 
to a value close to its initial value when x reaches 1.0. 
This is perhaps unsurprising since cat involves transport 
quantities S, p and pxy with opposite trends vs. x- By 
Eq. ^ we may flnd the curve of Uxy vs. T by adding 
the curves e n and S tan 6h (Fig. Ob) and dividing by p- 
The data in Figs. |2t andO) show that these terms are 
opposite in sign for x > 0.3. With increasing x, their mu- 
tual cancellation suppresses Uxy strongly. In particular, 
at the largest x (0.85 and 1.0), the cancellation is nearly 
complete and axy is very small below 100 K, i.e. the ob- 
served Cat is nearly entirely from the AHE of the backflow 
current. For 0.1 < a: < 0.3, S tamOn is negligible and cn 
largely reflects the behavior of axy- We exclude from our 
study the undoped compound CuCr2Se4 because eAf and 
p'xy were not resolved at low T. These trends emphasize 
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the importance of knowing all 4 transport quantities, in- 
stead of just ejv, to discuss J'^ meaningfully. 
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FIG. 3: (Panel a) Curves of ejv vs. T below 150 K in 5 
samples with doping 0.1 < a; < 1.0 showing nominal T-linear 
behavior at low T. The slopes vary non-monotonically with 
X. Panel (b) shows the Hall-current term StanOn measured 
in the same samples. For x > 0.3, StanOH is opposite in sign 
from ejv [the symbol key applies to both (a) and (b)]. Panel 
(c) shows the sharp change in the p-T profiles in the samples 
with X = 0.85 and 1.0. At low T, p at 0.85 is metallic, but at 
1.0, p reveals hopping between strongly localized states. 

Finally, the derived curves of Uxy vs. T are shown in 
Fig. 21 In contrast to the non-monotonic behavior of e^r 
vs. X, Uxy varies linearly with T as Uxy — b{x)T + c, 
where the slope b(x) now decreases monotonically as x 
increases from 0.1 to 1.0 (Fig. 0)). In all samples except 
X — 0.25, the parameter c - probably extrinsic in nature 
- is close to zero within our accuracy. The dependence 
of the parameter b{x) = [axy{T) — axy{0)]/T on these 
two quantities is of main interest. Figure compares 
how b{x) and ~ determined |^ from pxy above Tc - 
vary with x. Whereas, at small x, the decrease in b(x) 
seems to match that of n/j, b{x) falls much faster to zero 
at large x. 

A striking relation between them becomes apparent 
if we plot one against the other. Figure [3 shows that, 
above a threshold doping near 0.85, &(x) increases as a 
fractional power of nh — riho with nho a threshold density. 
This is consistent with b{x) increasing like the density of 
states (DOS), viz. b{x) - A/>. The DOS for the free- 
electron gas Afp (dashed curve) has a slightly stronger 
curvature than the data. Interestingly, the occurence of 
the threshold doping at x — 0.85 accounts well for a puz- 
zling change in the resistivity behavior when x exceeds 
0.85 (Fig. 13:). In general, slowly increasing x causes 
the resistivity profiles to change systematically, reflect- 
ing slight decreases in both and io (mean- free-path) . 
However, between 0.85 and 1.0, the change is sudden and 
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FIG. 4: [Panel (a)] Curves of a^y vs. T obtained by adding 
ejv and Standi/ (Eq. 0. The slope b{x) now falls monotoni- 
cally as X increases to 1.0. Panel (b) compares how the slope 
b{x) — Aaxy/T and nj, vary with x. 
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FIG. 5: Plot of b{x) — Aaxy/T against rik showing that 
above the threshold doping {x ~ 0.85), b{x) increases as a 
fractional power of Uh ~ nno- The dashed line is Aaxy/T — 
gAfp with g = 9.77 x 10"^" in SI units. 



striking. At a; = 0.85, p is T independent below 100 K 
consistent with a disordered metal. By contrast, at 1.0, p 
rises monotonically with decreasing T (Fig. |3t) . Between 
300 and 4.2 K, p increases from 6.3 to 32 mficm. At low 
T, conductivity proceeds by hopping between strongly 
localized states in an impurity band. Figure |5l confirms 
that we reach the extremum of the hole band near x = 
0.85. Further removal of carriers {x 1) affects states 
within the impurity band. 
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Knowing and p at each x, we may determine the 
mean-free-path £q in the impurity-scattering regime. Be- 
tween X = 0.1 and 1.0, decreases by a factor of 40. This 
steep decrease has no discernible influence on b{x). Com- 
bining these factors then, we have axy = gTAfp where 
g is independent of £o- We may boil down axy to the 
measurement of an "area" A by writing 

a., = ^^A/>, (r«Tc). (4) 

with kB Boltzmann's constant and e the electron charge. 
The value of g gives A = 33.8 if we assume Afp f^p- 
As the anomalous Hall heat current produced by E||x is 
= ayxTE, it shares the simple form in Eq. 0] The 
ratio Jy / Jy ^ T"^, as expected for a Fermi gas. 

We briefly sketch the anomalous- velocity theory 0, 0, 
1^. In a periodic lattice, the position operator for an elec- 
tron is the sum x = R + X(k) where R locates a unit 
cell, while X(k) locates the intracell position A fi- 
nite X(k) implies that x does not commute with itself. 
Instead, we have 0| [xj,Xk] = ie^'^™-Q,n with e-''^'" the 
antisymmetric tensor, which implies the uncertainty rela- 
tion AajjAa;^ ~ ri. The 'Berry curvature' r2(k) = VkxX 
is analogous to a magnetic field in k space [l^. In the 
presence of E, H adds a term that is the analog of the 
Lorentz force to the velocity Vk, viz. 

^Vk = Ve(k) -Ex fl{k). (5) 

The anomalous velocity term in Eq. [S] immediately im- 
plies the existence of a spontaneous Hall current 3h = 
— 2e^j^/£E X fifk) where is the unperturbed dis- 
tribution H, 0, H 0, 01 . The unconventional form of 
the current - notably the absence of any lifetime de- 
pendence - has made the KL theory controversial for 
decades 0, 0. However, strong support has been 
obtained from measurements of Lee et al. 8] show- 
ing that the normalized AHE conductivity cr'^y/uh in 
CuCr2Se4_xBrx is unchanged despite a 1,000-fold in- 
crease in p. 



In general, the off-diagonal term axy is related to the 
derivative of axy at the chemical potential p, viz. axy = 

^ ^ M ■ Using the result that a' is hnear 

in Hfi but independent of £q, and [drih/de]^ = Afp, we see 
that axy ~ TJVp, consistent with Eq. IH [By contrast, 
we note that the skew-scattering model 1^ would predict 
cr'xy ^ rihio and axy ~ TAfpio] 

Finally, A in Eq. |3| has the value 34 A . If Eq. El is 

indeed the origin of axy, A must be roughly the scale 
of O ~ AxjAxk- Hence the value A ^ ^ the unit-cell 
area seems reasonable (the lattice spacing here is 10.33 
A) . While a quantitative comparison requires knowledge 
of r2(k) over the Brillouin zone, the simple form of Eq. 0] 
seems to provide valuable insight on the anomalous heat 
current. 

A previous calculation of the Nernst coefficient was 
based on the "side-jump" model 0|. On scattering 
from an impurity, the carrier suffers a small sideways 
displacement S to give on average tan6'_f/ = S/£q. This 
was used to derive Qs ^ T(kp£)~^ . In our experiment, 
kp£ falls monotonically, with increasing x, while cn rises 
to a broad maximum near 0.25 before falling. Hence 
our experiment is in essential conflict with the side-jump 
model. From earlier experiments 41], an empirical form 
Qs = — (a + b'p)T has been inferred (a, b' are constants). 
This is not borne out in our data. 

Combining Nernst, Hall and thermopower experiments 
on the ferromagnet CuCr2Se4_xBrx, we have determined 
how axy (hence J*^) changes as a function of and r. 
At low T, we find that axy follows the strikingly sim- 
ple form axy ~ ATAfp, consistent with the anomalous- 
velocity theory for the AHE (Fig. In addition, a 
direct relation (Eq. ^ between M and axy is observed 
in the paramagnetic regime above Tq- 
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